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Abstract — Insulated gate bipolar transistors are the power 
semiconductors used for high current applications as a 
switching device. It is widely used in electrical and hybrid 
vehicles. It has become increasingly important to understand the 
reliability of these modules. The lifetime prediction is based on 
the assumption that the solder interconnections are the weakest 
part of the module assembly and that the failure cause is the 
inelastic deformation energy accumulated within the solder 
material. In this paper, the effects of thermo-mechanical fatigue 
of the solder layer interface have been investigated. A 2D model 
of 6 pack, 1200V IGBT module is used for the analysis. This 
study presents simulation of crack initiation and propagation 
under thermo-mechanical loading. Successive initiation 
technique is used to propagate crack in solder layer. A cyclic 
creep-fatigue damage under thermal loading is modeled using 
Energy Partitioning model to predict the fatigue life of an IGBT 
module. 

Keywords — Successive initiation technique, IGBT, Energy 
partitioning model, Finite element analysis. 


Battery Management: The battery management system 
controls battery state during charging and discharging. 
Intelligent functionality is needed to extend the battery 
lifetime, which has a considerable impact on the total cost of 
ownership. The State of Health (SoH), State of Charge (SoC) 
and Depth of Discharge (DoD) of the battery is permanently 
monitored. 
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i. Introduction 

The Now a day’s many companies are working on hybrid 
vehicle and electrical vehicle. Main Inverter, DC/DC 
Converter, Auxiliary Inverters/Converters, Battery 
Management and On-board Charger are the components of 
HEV and EV. Figure 1 shows all the components. 

Main Inverter: With an electric drivetrain, the inverter 
controls the electric motor. This is a key component in the car 
as, similar to the Engine Management System (EMS) of 
combustion vehicles, it determines driving behavior. Not only 
does the inverter drive the electric motor, it also captures 
energy released through regenerative breaking and feeds this 
back to the battery. As a result, the range of the vehicle is 
directly related to the efficiency of the main inverter. 

DC/DC Converter: Different voltage levels are required by 
the various electronic components in an EV. High-voltage 
batteries with different voltage levels are currently available 
on the market. In addition, the power classes scale from lkW 
to 5kW depending on the number of low- voltage applications. 
In EVs and HEVs, the DC/DC converter supplies the 12V 
power system from the high voltage battery. 

Auxiliary Inverters/Converters: Power on demand 
increases the efficiency of electric vehicles. With HEVs and 
EVs, former belt-driven devices classified as auxiliary drives 
have been electrified and integrated into the power system to 
deliver power on demand. Typical auxiliary systems supplied 
from the high-voltage battery include air conditioning, 
electronic power excellent thermo-mechanical stability with 
effective thermal management of IGBT chips. 


Figure 1 . Shows components of HEV 

On-board Charger: The battery in an electric vehicle is 
useless without a battery charger. And all electronic systems 
depend on the battery for power. With an on-board charger 
unit, the battery can be charged from a standard power outlet. 
This is very important function of electrical vehicle 

Power modules contain one or more semiconductor 
switches packaged together for easy connectivity. Power 
semiconductor modules play key role in power electronic 
system. The module can be used to enclose thyristors, diodes, 
FET, IGBT, or similar semiconductor switches. Nowadays, 
IGBTs are widely used as a result of good switching 
performance combined with fairly low conduction losses[l]. 
IGBT is a four layer semiconductor device. It combines the 
voltage characteristics of a bipolar transistor and the drive 
characteristics of a MOSFET. The popularity of the IGBT has 
increased in recent years due to an increase in high voltage, 
high power application. They are available in the range of 
voltage rating from 300 to 1200 volts and current rating from 
15 to 100 amps for a single die. IGBT modules have current 
rating well into the 100’s of amps [2]. 

IGBT power modules are used in applications like 
automotive, traction, solar inverters. In hybrid electric 
vehicles, an electric power conversion system that includes an 
inverter and converter. They are used to convert the power 
generated by the engine into electrical energy, to charge and 
discharge a battery, and to drive the motor. The electric power 
conversion system typically uses an IGBT module as its main 
switching device. The failure modes that have been identified 
are A1 wire bond failure and crack in the solder layer between 
the base plate and substrate. Hence, the major concern from 
the reliability standpoint of view is the effect of thermal 
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cycling on the insulating substrate to the base-plate solder 
interface. In this study, fatigue life prediction for solder 
interconnects in IGBT module by using E-P damage model 
and the SI method is done. 

Ladani et al have discussed the use of SI method for 
lead-free BGA solder interconnects [3]. In this research an 
experiment was designed and implemented on printing and 
reflow processes of Pb-free solder interconnects in order to 
obtain error-seeded specimens. The specimens were 
characterized for voids and subjected to temperature cycling 
in order to evaluate the durability of specimens under 
thermo-mechanical loading. Statistical analysis of the void 
characterization as well as durability suggested the 
detrimental effect of voids. An existing E-P model was 
modified by SI method and used to explicitly predict damage 
initiation and propagation in BGA solder joints. Two 
techniques were used to predict the life, averaging technique 
and successive initiation [3]. James Gyllenskog has reported 
the use of E-P model coupled with a SI technique to see the 
initiation of crack and fatigue life prediction of the T-38 
aileron lever for materials aluminum 2014-T6 and 7050-T74 
[4]. 

E. BACKGROUND 

Transistor is a three lead semiconductor device that acts as an 
electrically controlled switch or a current amplifier. BJT is a 
current controlled device, whereas FET and MOSFET are 
voltage controlled devices. The simple addition of an extra 
P-N junction to the drain of the MOSFET changes the 
unipolar device into bipolar junction transistor. Combination 
of an insulated gate input and bipolar output makes the IGBT 
an excellent power switch for medium frequency and high 
voltage application. It combines features of both of these 
devices. Power electronics packaging technology has been 
developed for several generations, involving material 
upgrading, structure improvement and interconnection 
technique innovation. IGBT is presently one of the most 
popular device due to its wide ratings and switching speed of 
about 100 KHz, which makes it an easy voltage drive. Among 
the new power devices, IGBT devices are being more 
accepted and increasingly used in traction application such as 
locomotive, elevator, tram and subway. 

Power semiconductor modules play key role in power 
electronic system. The module can be used to enclose 
thyristors, diodes, FET, IGBT, or similar semiconductor 
switches. The reliability of solder interconnects has always 
been an area of concern for microelectronics applications and 
continues to be a concern for power electronics as well. The 
study has been done by Hua Ye et al[5], on different failure 
modes in power electronic devices. They concluded that 
thermal cycling has significant effect on the solder layer 
reliability. They also stated that using finite element analysis 
and damage mechanics constitutive model it is possible to 
predict number of cycles to failure for solder layers. Sun et 
al[6] proposed CZM-based thermo-mechanical fatigue 
model. This model was developed to simulate the 
delamination process at the die-solder interface and to predict 
the lifetime of solder joint in IGBT module. In this model the 
path of the crack is known. To analyze this crack propagation 
process by delamination, a very thin layer of the solder joint at 
the die-solder interface is replaced by a layer of cohesive 
elements with the same thickness. The damage will be 


concentrated on this cohesive layer which allows simulating 
the progressive crack initiation and propagation taking place 
at the die-solder interface. The thickness of the cohesive layer 
will influence more or less the thermo-mechanically induced 
strain and stress fields. Hence, the effect of the thickness on 
the rate of damage evolution and the lifetime is studied. There 
are several different approaches like stress based approach, 
strain based approach, energy based approach and E-P 
damage model approach that are used to predict the life of 
solder joints. Present study focuses on fatigue life prediction 
for solder interconnects in IGBT module by using SI 
technique and E-P model. SI technique is used to propagate 
crack in solder layer interface between substrate and base 
plate. E-P damage model approach is used to predict the 
lifetime of an IGBT module. E-P damage model assumes that 
cyclic fatigue damage is due to combination of elastic, plastic 
and creep deformation mechanism [3]. 

m. MODELING METHDOLOGY 

A 2D IGBT module is used for the analysis. Module 
contains die soldered onto a DBC substrate which is made out 
of Aluminum Nitride (AIN) material for internal electrical 
insulation. It is then soldered to the base plate. Figure 1 shows 
a 2D view of the multi-layer IGBT structure. Dimensions of 
each layer within the module construction are given in 


Table 1 






copper 2 


substrate 

copper 1 

solder layer 1 

base plate 


Figure 1: Schematic diagram of a conventional IGBT 
module [7]. 


Part 

Material 

Total Dimensions (mm) 

Die 

Silicon 

12.7 x 12.7x0.42 

Solder Layer 2 

SAC 305 

12.7 x 12.7x0.08 

Copper 2 

Copper 

41.7x27.5x0.3 

Substrate 

AIN 

41.7x27.5x0.6 

Copper 1 

Copper 

41.7x27.5x0.3 

Solder Layer 1 

SAC 305 

41.7x27.5x0.21 

Base Plate 

AlSiC 

174.1 x 92x3 


Table 1. Dimensions of IGBT module 
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Part 

Material 

Young 

Modulus 

(Gpa) 

Poisson's 

ratio 

CTE 

(ppm/°c) 

Die 

Silicon 

131 

0.3 

2.8 

Solder 
Layer 2 

SAC305 

MISO Properties 

20 

Copper 2 

Copper 

121 

0.3 

17.3 

Substrate 

AIN 

330 

0.24 

4.5 

Copper 1 

Copper 

121 

0.3 

17.3 

Solder 
Layer 1 

SAC305 

MISO Properties 

20 

Base Plate 

AlSiC 

188 

0.29 

8.75 


analysis. To apply boundary condition, all nodes at bottom 
surface of base-plate are fixed. 

Below shows thermal time history for -55°C to 160°C 
temperature range. The results are calculated at every thermal 
cycle and actual results consider are difference between last 
thermal and second thermal cycle. This technique gives us 
accurate results. Temperature cycle having time period 15 
minutes ramp down and ramp up, 15 minutes dwell. Total 
thermal cycle required 60 minutes to complete. Same thermal 
cycle consider for temperature range 0°C to 100°C 


Table 2. Material properties used for IGBT module 

Mesh sensitivity analysis was carried out where plastic 
work varies in the range of +/- 2%. Figure 2 shows 2D finite 
element model of IGBT module. Element plane 183 is used 
for the analysis. 




Base plate 


Figure 2: Two dimensional finite element meshed structure 

Material properties used for the analysis are given in Table 
2. The solder layer is modeled using multi-linear isotropic 
properties and creep properties. The primary creep is 
neglected due to negligible damage in thermal cycling. 
Secondary creep is considered for the analysis. For solder 
materials Garofalo constants are used which are provided in 
Table 3 [8]. 


Figure 3: Thermal time history. 

IV. ENERGY PARTITIONING DAMAGE MODEL 

A cyclic creep-fatigue damage under thermal cycling is 
Tnodeled using an E-P model. This model predicts cyclic 
creep fatigue damage based on deviatoric energy densities U e , 
W p and W c for a typical load cycle. The damage due to each of 
these deformation mechanisms is determined by using a 
power law as provided in equation (1) to (3). 

, U e =U e0 N e b .... (1) 

1 W p =W p0 N p c .... (2) 

W c =W c0 N c d .... (3) 

Figure 4. shows log-log plot of elastic, plastic and creep 
energy density (i.e. U e o, W p0 , and W c0 respectively ) versus 
cycles to failure. The exponents b, c, and d are their 
corresponding slopes. W p and W c are obtained from FEA 
analysis. Energy densities are calculated on those elements 
which are eliminated in every run. 


Solder alloy 

Creep properties 

SAC305 

a (MPa)' 1 

n c 

A 

Q(J-mol-K) 

0.18 

2.3 

6.07 

55793 


Table 3. Secondary creep Garofalo constants [8] 

IGBT module is subjected to thermal cycles ranging from 
-55°C to 160°C and 0°C to 100°C with 15 min ramps and 15 
min dwells. Figure 3 shows the thermal cycle range from 
-55°C to 160°C [5] . Plain strain condition is considered for the 



Figure 4: Work density verses number of cycles to failure 


[ 3 ] 
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Damage caused by each type of deformation can be calculated 
by substituting equations (1), (2) or (3) in equation (4). 
Damage is defined as reciprocal to number of cycles, as 
follows: 


1 



The total damage can be calculated by adding damage caused 
by elastic, plastic and creep deformation as seen in equation 
(5). Usually damage caused by elastic deformation is 
negligible. 

®tatal = -b Ai + D c 


1 

N f 




....(6) 

The total number of cycles to failure N f is calculated from 
equation (6) [3]. 


V. SUCCESSIVE INITIATION METHOD 

The SI method was first introduced by Okura[3]. SI 
analysis involves several steps for damage initiation and 
propagation shown in Figure 5. SI method is implemented 
with the help of finite elements. The damage initiation site is 
first identified with the help of an E-P damage model, for a 
typical cyclic loading condition. 





Inelastic energy includes plastic and creep work per cycle. 
This energy is calculated using FEA for all the elements to 
estimate the amount of damage by an E-P damage model. The 
elements with maximum damage can be identified from 
contour plots. The initiation zone is recognized by selecting 
all elements which exceed a suitable damage threshold. In 
order to minimize number of runs and to facilitate reasonable 
computational time, damage threshold in our case was 
selected to be 60% of the maximum value of damage. The 
propagation path for maximum damage evolution is then 
monitored. Elements in the damage initiation zone are 
removed from the structure by eliminating or killing the 
damaged elements. The cyclic loading is repeated on the 
remaining surviving structure and new damage zones are 
identified based on 60% of the maximum value of damage. 
The process is incrementally repeated, until a complete failure 
path is established and the critical solder joint loses its ability 
to carry any further loads [9]. 

The total number of cycles to failure is calculated using 
MATLAB algorithm for thermal cycles ranging from -55°C 
to 160°C and 0°C to 100°C. Calculations and tables shown 
below indicate values for thermal cycle -55°C to 160°C. 
Similarly, values have been calculated for thermal cycle 0°C 
to 100°C. Following is the procedure to calculate total 
number of cycles to failure. The term "Run" which is used in 
Table 4 to Table 10, refers to the thermal cycling applied on 
FEA model. Calculated damage is accumulated in individual 
elements. The threshold damage criteria is applied to kill the 
elements having value excess of defined threshold. Above run 
procedure is repeated till the complete delamination of solder 
layer is achieved. 

Calculated creep work values for the runs are the per cycle 
work (i.e. creep work developed at the end of last thermal 
cycles). Table 4 shows the values from run 1 to run 56. 
Calculated plastic work values for the runs are the per cycle 
work (i.e. plastic work developed in the last thermal cycle). 
Table 5 shows the values from run 1 to run 56. In following 
tables "0" value shows set of element having energy more 
than threshold. From values of creep work and plastic work 
estimate number of cycles due to creep and plastic work (N c & 
N p ) by using work density verses number of cycles to failure 
graph (Figure 4). 

Following equations are used to calculate damage due to 
creep and plastic work. 


1 



( 8 ) 

We get damage per cycle by addition of damage due to 
creep and plastic. Table 6 shows damage per cycle values. 


Figure 5: Damage propagation algorithm using SI 
method[9]. ^ 
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Reciprocal of damage per cycle gives value of discrete N f . 

Table 7 shows discrete N f values. 

1 


Discrete l V/ = 


Cycle 


(10) 


Second row of total damage is calculated by following 
equation. 

T otal Damage = 

[( Cumulative Nf — Cumu latim 

1 — I ■ J g Currs-n t row 


L 


X Damage per Cycle ^ 


illative Nf 

S 2/ I 

"it rove vahie 


Cumulative N f and total damage values are calculated 
simultaneously. First row of cumulative N f is equal to the first 
row of discrete N f . First row of total damage is equal to first 
column value of cumulative N f divided by current column 
value of cumulative N f . 

Second row of cumulative N f is calculated by following 
equation. 

(.1 — (X o tal D amag e) p rmv J 

Cumulative N f = ■ 


/Damage x \ 
V 'Cycle) 


current row 


,(ii) 


( Cumulative Nf ) 




Cl -0.704992] 


0.005516 


.( 12 ) 


Above equation (11) is used to calculate remaining 
cumulative N f values. Table 8 shows cumulative N f values. 
For e.g. 


Above equation (12) is used to calculate remaining total 
damage values. Table 9 shows total damage values. 

For e.g. 

Total Damage @ R[lY1 ^ 

= 1 - [(109.4635 - 53.4759) x 0.00370] 

= 0.792355 

Total Nfis to addition of diagonal values of cumulative 
N f [10]. Table 10 shows total N f values. 

For e.g. 


Total Nf 


pRun 2 


= 230.7731 4- 53.4759 


= 53.4759 


= 292.2491 


Table 4.Creep Work 



Table 5. Plastic Work 
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Table 6. Damage per cycle 



Table 7. Discrete N f 



Table 8. Cumulative N f 



Table 9. Total Damage 
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VI. RESULTS 

The fatigue lifetime is often represented by the number of 
cycles to failure. Fatigue damage in solder layer was found 
due to accelerated thermal cycle. The numbers of cycles to 
failure of the IGBT module are around 2266 cycles and 2400 
cycles for thermal cycles -55°C to 160°C and 0°C to 100°C 
respectively. Figure 6 shows the 2D IGBT module before 
crack initiation and Figure 7 shows module after crack 
initiation and propagation in solder layer for thermal cycle 
-55°C to 160°C. 



Figure 6: 2D IGBT module before crack initiation. 



Figure 7: 2D IGBT module after crack initiation and 
propagation. 


vn. CONCLUSION 

In this paper, 2D finite element analysis of a multilayered 
IGBT module under thermal cyclic loading is carried out. 
Damage is calculated in individual elements using E-P model 
coupled with SI technique. Using threshold damage criteria, 
the elements are killed. The crack is initiated at the outer side 
of the bottom solder layer interface in between base plate and 
substrate due to CTE mismatch. The crack is propagated 
towards the inner side of the bottom solder layer till it is 
completely delaminated. As a result delamination of solder 
layer, IGBT module fails to carry any further mechanical and 


electrical loads. Parametric study performed for IGBT 
module. Considered thermal cycles -55°C to 160°C and 0°C to 
100°C, we get number of cycles to failure are 2266 and 2600 
respectively. Thus fatigue life is predicted for solder 
interconnects in IGBT module. Conclusion of above 
parametric study is thermal cycle AT value decreases the 
number of cycles to failure increases. Future work will 
involve fatigue life prediction of the IGBT module by 
comparing various solder alloys with different thermal cycles. 
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